Raft Instability of Biopolymer GelsS 
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Following recent X-ray diffraction experiments by Wong, Li, and Safinya on biopolymer gels, 
we apply Onsager excluded volume theory to a nematic mixture of rigid rods and strong "7r/2" 
cross-linkers obtaining a long-ranged, highly anisotropic depletion attraction between the linkers. 
This attraction leads to breakdown of the percolation theory for this class of gels, to breakdown of 
Onsager's second-order virial method, and to formation of heterogeneities in the form of raft-like 
ribbons. 



87.15.-v (Biomolecules: structure and physical 
models and theories of liquid crystal structure), 

The demonstration by Lars Onsager in 1949 M that 
solutions of long thin rods undergo a first-order phase- 
transition from an isotropic to a birefringent nematic 
phase has remained a landmark achievement of the- 
oretical statistical physics. He showed that at the 
isotropic/nematic transition point the rod volume frac- 
tion <j> is surprisingly low, of the order of the aspect ratio 
D/L of the rods, D being the rod diameter and L the rod 
length. Stiff biopolymers provide interesting applications 
of Onsager theory. In particular, the rod-like filamentous 
protein actin, which carries many biophysical functions 
PL exhibits an isotropic-to- nematic transition at a 
critical volume fraction in approximate agreement with 
the Onsager criterion ||. Inside cells, micron-size actin 
filaments are part of a gel structure, the cytoskeleton, 
with the filaments cross-linked by reversibly bound pro- 
teins like a-actinin. The unusual elastic properties of 
in-vitro actin gels j|,D have recently been the focus of an 
intense theoretical effort JjJ . 

In-vitro studies of the sol-gel transition of fixed-length 
actin filaments are currently interpreted in terms of per- 
colation theory [||, in which linkers are assigned, at ran- 
dom, to rods providing connections to neighboring rods 
with no preferred crossing angle. The sol-gel point, where 
the real and imaginary parts of the elastic moduli are 
comparable in magnitude, is identified as the percolation 
threshold, i.e., the linker concentration for which a con- 
nected "path" of linked rods stretches across the sample 
for the first time. An important validity condition is that 
the concentration of "native" rod-rod contacts - i.e., the 
mean number density of rod-rod contacts before linkers 
are added - exceeds the percolation threshold, otherwise 
the reduction in configurational freedom of the rods im- 
posed by the linkers would lead to phase-separation. The 
number of native contacts per rod in a dilute isotropic so- 
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lution of rods is of the order of (L/D)d> so a rod-linker 
solution with <\> of order D/L is expected to exhibit a sol- 
gel transition that is reasonably well described by per- 
colation theory. The same argument indicates that the 
structure factor 5(q) should resemble that of a pure On- 
sager nematic. 



(a) Nematic 




(b) n/2 Gel 




Charged Rods 



(c) A single raft (d) A stack of rafts 

FIG. 1. Mixtures of rigid rods in the nematic phase and "ir/2" 
linkers display a complex range of structures, (a) At very low linker 
concentrations and binding energy, excluded volume effects prevent 
cross-linking, (b) Second-virial theory predicts "7r/2" tetratic or 
cubatic gels at high linker concentrations, (c) At low linker con- 
centrations a long-range depletion attraction leads to heterogeneous 
structures in the form of raft-like ribbons along the optical axis, (d) 
Rafts can stack on top of each other as a result of inter-raft attrac- 
tion leading to phase separation of the linker-rich stacks from the 
nematic phase. 
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Recent high-precision low-angle synchrotron X-ray 
studies § of the microscopic properties of mixtures of 
actin filaments and positive divalent counterions report 
results that are in conflict with these expectations. When 
divalent counterions are added to a nematic phase of 
actin filaments, a birefringent gel phase is observed at 
linker concentrations that depends sensitively on the (av- 
erage) rod length L, whereas percolation theory would 
predict that, beyond a certain length, the critical linker 
concentration is determined only by the rod concen- 
tration. In addition, a harmonic sequence is found in 
the structure factor 5(q) of the gel phase in the long- 
wavelength range 1/L < q < 1/D along directions ap- 
proximately perpendicular to the optical axis. This is 
incompatible with the Onsager nematic [[To) , but it is typ- 
ical of layered structures, such as smectic liquid crystals. 
The X-ray experiments on actin gels as well as similar 
results obtained on other biopolymers such as DNA Jll]] 
clearly conflict with a simple "Onsager-Percolation" type 
description. 



It is the claim of this paper, that the unusual features 
of biopolymer gels can be physically interpreted if we al- 
low for the fact that linked biopolymers have preferential 
crossing angles. First, actin linker proteins are known 
to impose different crossing angles, depending on their 
molecular structure. Second, biopolymers like actin or 
DNA carry a high linear charge density, which allows 
them to be water soluble. The electrostatic repulsion 
between two uniformly charged rods linked at a single 
point depends on the relative angle 7 as f/| sin7| [^2[ , so 
the optimal crossing angle between two linked biopoly- 
mers is, in general, expected to be large. Based on this 
consideration, we studied - as a simple model for biopoly- 
mer gel formation - a dilute mixture of rigid rods with a 
concentration p v of the order of 1/DL 2 , in the presence 
of a concentration pi of freely sliding, reversible linkers 
connecting pairs of rods at 7r/2 crossing angles (binding 
energy eq). When pi = 0, second- virial theory predicts 
a first-order phase transition from an isotropic to a ne- 



matic phase. It is easy to demonstrate |13 that at non- 
zero linker concentrations, second- virial theory predicts a 
first-order phase transition (for pi oc p p ) from the nematic 
phase to a birefringent "7r/2" gel phase with tetragonal 
symmetry Q (see Fig. 0b). 




FIG. 2. Schematic view of (a) the depletion volume of a trans- 
verse rod in an Onsager nematic and (b) the overlap of two deple- 
tion volumes surrounding two transverse rods for R. = (x, 0, 2) and 
Yti = 0. 

To examine gel formation beyond second-virial the- 
ory, we follow what happens when linkers are added, one 
by one, to the rod solution. A single linker connect- 
ing two rods produces a cross-like structure, with the 
two rods free to slide with respect to each other. It is 
easy to show that the free energy has a minimum when 
one arm of the cross extends along the optical axis (the 
z-direction) with the other arm taken to be the trans- 
verse x-direction [lq ]. As shown in Fig. ^a, the trans- 
verse rod creates an anisotropic depletion volume 8V\ 
for the centers of mass (CM's) of the free rods. This 
depletion volume depends on the angle 7/t between the 
axis of a free rod and the axis of the transverse rod by 
the Onsager result SVi(-ff t ) — 2DL 2 \ sin7/t|. The en- 
tropic self-energy of the linked cross is computed as the 
osmotic work 5W\ ~ Tl osm (8Vi('jft)}g required to remove 
the CM's of the free rods from the depletion volume, with 
Ilosm the osmotic pressure of the nematic and with ( )g 
denoting an orientational average. The effect of the arm 
of the cross along the optical axis is smaller by a factor 
D/L and will be neglected. Using Onsager theory for the 
orientational distribution of the free rods, we find: 



SWi 



2 C/ (1 4 
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ci) 

- A6> 2 /4) 



(Isotropic, ci < 3.3) 
(Nematic, c^ > 4.5) 



(1) 



with A8 2 being the (small) variance of the polar angle 
with respect to the optical axis, and c = (ir/4)DL 2 p p 
the Onsager dimensionless rod concentration. At the 
isotropic/nematic transition point, the dimensionless free 
energy cost SWi/ksT ~ 33.2 assumes a large value inde- 
pendent of the dimensions of the rod, or the rod concen- 
tration. 

Now add a second linker with a second rod making a 
7r/2 angle with the optical axis (see example in Fig. ||b). 
Let the relative distance between the CM's of the two 
transverse rods be R, let the angle of the second rod 
with the x direction (i.e., the relative angle) be jtt-> an d 
let <5V2(R, 7tt) be the excluded volume that is shared be- 
tween the two rods. The entropic free energy cost of the 
two-rod system will then be 



<Wa(R,7«) - 2<Wi - U osm (SV 2 (R^tt))e 



(2) 



In the limit R = |R| — > and jtt = 0, the two rods coin- 
cide so 5W 2 (0, 0) ~ SWi. On the other hand, for R > L, 
there is no shared excluded volume so <5Wa(R, 7«) — 
2SWi. It follows that the two linkers attract each other 
with a potential energy that has a typical range of or- 
der L and a typical binding energy of order SWi. This 
effect is related to the well-known depletion attraction 
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Jig , the attraction between large objects in a surround- 
ing medium of small objects driven by the increase in 
configurational entropy as the larger objects coagulate. 
The difference is that in the present case an effective 
attraction between small objects (linkers) is generated 
by an increase in configurational entropy of large objects 
(rods). Using once again the Onsager orientational distri- 
bution, we find that, according to Eq. ||, the linker-linker 
depletion attraction drops off inversely proportional to 
the distance between the rods: 



SW 2 {z, ltt )-2SW 1 
6W1 



D 



-at 



A6\z\ 



D 



L\ smj tt \ 



1 



(3) 



where z is the separation between the rods along the 
optical axis and a\ ~ 6/yfn. The interaction is highly 
anisotropic: for small \z\, there is a strong torque be- 
tween transverse rods favoring parallel alignment (i.e., 
Itt = 0). Moreover, the mechanism operates only as 
long as the CM of the second rod lies below or above 
the first rod within the two narrow wedges bordered by 
\y\ ~ A9\z\ centered on the first rod. 

Equation § is valid for D <C A9\z\ <C L. For small 
separations A0\z\ <C D the potential has a linear depen- 
dence on z: 



5W 2 {z,0) - 26Wi 



SWi 

- U-02- 



A6\z\ 
~D~ 



D 



L\ sin7 M 



(4) 



with a,2 — f/2 v / 7r. This result can be understood by 
noting that for perfectly aligned rods the shared deple- 
tion volume of two transverse rods would have a linear 
dependence on separation JT?]] . 



(a) 




8n v <0 



8n„>0 



(b) 




8» >0 



hn v <0 



FIG. 3. (a) Schematic view of the long-range depletion inter- 
action between two transverse rods. The continuous lines indicate 
the director field, (b) The quadrupole symmetry of the distortion 
zone around a transverse rod. 

In the opposite limit of rod separations large compared 
to the rod length L, the continuum description of a ne- 
matic Jl8[ can be used in terms of a locally varying direc- 
tor field n(r) (see Fig. ||a). The perturbation of n(r) far 



from the fixed rod can be written as n(r) ~ (0, 6n y (r), I), 
where Sn y <C 1. The leading term in the elastic free en- 
ergy in this geometry is the splay term: 



splay 



§i J d 3 r(V5 % ) 



(5) 



where Ku is the splay elastic constant and is equal 
to K n ~ 0.06L»L 4 p2/c B T for rigid rods @. Treating 
Sn y (r) as a variational function and minimizing Eq. § 
leads to the Laplace equation V 2 5n y = 0. As demon- 
strated in Fig. ^3, the appropriate solution of the Laplace 
equation must have the symmetry of a two-dimensional 
quadrupole so Sn y (y,z) = Qyz/(y 2 + z' 2 )' 2 . The value of 
the quadrupole moment Q is determined by the condition 
that the typical deviation Q/L 2 of Sn y (r) for y 2 + z 2 ~ L 2 
should match the angular deviation D/L of the free 
rods as obtained from excluded volume arguments (see 
Fig. ||a). This condition leads to Q cx DL times a func- 
tion of AO. From the electrostatic interaction energy be- 
tween quadrupoles in two dimensions, Q 2 cos(4c/?)/r 4 (in 
polar coordinates r, ip) , it follows that the interaction en- 
ergy per unit length between two parallel (infinite) rods 
must have the form: 



^splay 

Lk R T 



g(A9) cos(4t^)c 2 



,DL 2 



r > L 



(6) 



where the function g(A9) includes the dependence on the 
nematic order parameter. The deformation energy F sp \ a y 
should be matched with the right hand side of Eq. g 
at distances r ~ L yielding g cx I/c. It follows from 
Eq. ^ that the anisotropy of the linker-linker interaction 
is significantly less pronounced in the regime r L. 
Moreover, the interaction is repulsive when the separa- 
tion vector between the CM's of the rods makes an angle 
near 7r/4 with the optical axis. The various regimes are 
shown schematically in Fig. 

bW/bW^ -2 




-1 



intermediate 

(~l/z) 



continuum 
limit 



(~1/Z 4 ) 



FIG. 4. The depletion attraction between the two rods for 
y = 0. 
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We now turn to the implication of these results. Ac- 
cording to Eq. [|, the characteristic energy scale 8W\ of 
the long-range depletion attraction between 7r/2 linkers 
is large compared to the thermal energy ksT. This im- 
plies that percolation theory is not applicable since link- 
ers will cluster even at very low concentrations. Next, 
second-virial theory is not reliable cither since it gives 
no indication of this form of phase separation. Indeed, 
the cross-links resemble flat square discs of dimensions 
L x L and second-virial expansions are known not be 
accurate for mixtures of rods and discs [^0|. The large 
value of SWi/ksT would have another important con- 
sequence: unless the linker binding energy is sufficiently 
high, crosses simply would not be able to form. Assuming 
chemical equilibrium between linkers in solution, linkers 
adsorbed on isolated chains, and linkers forming a cross- 
structure, we find that the fraction / of linkers able to 
form a cross is of order: 



/ 



l + (Tr/2c)e^o-sw 1 )/k B T 



(7) 



This means that cross-links disappear when 5W\ is 
large compared to £$, i.e., when rods are so long that 
DL 2 3> eo/ksTpp. This is consistent with the experi- 
mental observation that the critical linker concentration 
for actin gelation rapidly increases with rod length. For 
an actin solution with p p ~0.1mM, a linker concentra- 
tion pi ~lmM and Eq ~ 15ksT, the rod length should 
not exceed a value of order lOOnm. 

What would be the structure of the clusters? Accord- 
ing to Eqs. |3[ |], the depletion attraction between two 
transverse rods is highly anisotropic at short distances 
and maximal if the two rods are parallel and located just 
below each other with a common plane whose normal 
is perpendicular to the optical axis. Adding additional 
linkers and transverse rods in the same way produces an 
elongated ribbon structure resembling a disordered raft 
(see Fig. |l|c) . The sol-gel transition point should be char- 
acterized, not by percolation, but by the formation of 
elongated, ribbon-like rafts extending through the sam- 
ple along the optical axis. Because two parallel plates 
immersed in a nematic attract each other |2lf| - with a 
force per unit area that depends on the inter-layer spac- 
ing h as ksT/h 3 for large h - depletion attraction should 
also operate between different rafts, so the ribbons could 
be several layers thick (Fig. |l]d) . A multilayer ribbon-raft 
indeed would have along the transverse direction a struc- 
ture factor 5(q) similar to that of of a smectic, although 
without long-range order. Along the optical axis, there 
would have to be a peak in S(q z ) when 2n/q z equals the 
spacing between transverse rods inside the raft, which is 
indeed present jnj . 

Our method only applies to a small number of link- 
ers. At finite linker concentration, rafts and ribbons are 
expected to act as long-lived kinetic intermediates. The 
actual thermodynamic equilibrium state is expected to 



be either a ir/2 gel or a biaxial nematic JlO| , pl| , in phase- 
coexistence with linker-poor solution. Numerical simula- 
tion methods of the model should be able to address this 
issue. Next, we did not allow for the fact that linkers may 
prefer the endpoints of rods nor deviations away from 
ir/2 for the preferential crossing angle. Finally, actin cy- 
toskcletons do appear to resemble tt/2 gels p2| with no 
evidence of raft formation. Actin filaments inside cells 
do not have a fixed length however: active polymeriza- 
tion and depolymerization processes are constantly tak- 
ing place and in-vitro experiments of gel formation by 
such "living" actin filaments would be of great interest. 
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